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THE DETERMINATION OP SPAN LOAD DISTRIBUTION AT 
HEGH SPEEDS BST USE HIGH-SPEED WIND- 
TUNNEL SECTION DATA 
By John Boshar 



SUMMARY 



A tabular method la presented for determining the 
span load distribution at high Mach numbers by utilizing 
high-speed airfoil section data. The method, designated 
the generalized method. Is an easily applied process of 
successive approximations by which a general application 
of the llftlng-llne theory may be used to determine the 
span load distribution for wings composed of sections 
having arbitrary lift curves. Aii example Is given to 
show how this method Is used. A comparleon of span 
load distribution obtained by the generalized method 
using high-speed data is made with results obtained by 
the strip-theory method using high-speed deta and by 
the conventional method of applying liftlng-ilne theory, 
which utilizes low-speed data. 

The results of the computations Indicate that the 
loading changes associated with Mach nvimbcr may be 
great enough to require modification in the current 
method of computing span loading for design purposes. 



INTRODUCTION 



in some recent high-speed airplane flights a niimber 
of accidents and near accidents have occurred, which 
coxild be associated, in part, with changes in wing span 
load distribution at high Mach numbers. These changes 
have been manifested by both inordinate changes in air- 
plane stability and formation of permanent wing wrinkles 
at load factors lower than should be expected on the 
basis of static-test results. 

Qualitative considerations of the effect of Mach 
nximber on two-dimensional airfoils indicate that some 
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changes In loading are to be expected In the usual 
operating lift range, because the thicker Inboard sec- 
tions of a typical wing vould experience a con^X'easl- 
blllty stall earlier than the outboard sections and would 
consequently require the outboard sections to carry a 
greater part of the load If the lift is to be maintained. 
Except for some unpublished high-speed v;lnd-tunnel re- 
sults of wake maasureraents behind a tapered wing, which 
verified the conclusion that stalling occurs earlier 
on the root sections than on the tip sections, no direct 
experimental data exist on the subject and the actual 
magnlt\xde of the span-loading changes has been question- 
able. 

Although designers are aware, therefore, that 
changes occur In the soction 11 ft curves after the occur- 
rence of a corcpressiblllty stall, airplane wings are 
still built to carry the limit loads distributed in ac- 
cordance with a llftlng-llne theory that Includes the 
assumption that the individual sections alon^ the span 
have a constant llft-cjirve slope throughout the entire 
operating range. 

The purpose of the present report Is: 

(1) To present a tabular method oy. which the 
lirtlng-llne theoi'y may be easily applied to the deter- 
mination of the span load dleirioi;tlon of a wing, re- 
gardless of the type of section lift curves. 

(2) To present, for two hypothetical wings, com- 
parisons of the span load distribution as determined by 
means of 

(a) The strip theory, in which high-speed 

wind-tunnel data are used 

(b) The conventional application of the 

llftlng-llne theory, in which low- 
speed wind-tunnel data are used 

(c) The generalized method of applying lifting- 

line theory, in which high-speed" wind- 
tunnel data are used 
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(3) present an example showing how the general- 
ized method may be used to determine span load distribu- 
tion at high speeds. 

S7MB0LS 

A aspect ratio (b^/S) 
b wing span 

section lift coefficient 
Cj^ wing lift coefficient 
^BM bendlng-moment coefficient 
c chord at any span 
o* mean chord 
M Vach number 

m, n, r, s constants deflnl^ig limits of summations 

n load factor 

q dynamic pressure 

S wing area 

V airplane speed 

w downward component of velocity 

W airplane weight 

y distance from wing root, semlspans 

a angle of attack; with subscript g, geometric; with 
subscript e, effective 

e downwash angle 

^ constant coefficients In downwash equations 

A prime used with a symbol Indicates a partlqular 
value of the quantity. 
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MEI'HODS OP COMPUTING SPAN LOAD DISTRIBUTION 



A number of methods have been used to determine 
spanwlae distribution of air load on airplane wings. 
The most prlmj.tlve of these, the strip method, has been ■ 
used to obtain rough solutions but. Inasmuch as this ... 
method neglects the presence of downwash. It has not 
been consjdered accurate enough for design purposes. 
The method that Is used In desl^ accounts for the 
effect of dovmwash by an application of the liftlng-llne 
theory.. This method which, for convenience. Is referred 
to herein as the "conventional metliod," has been made 
adaptable to direct computation for the special case of 
lift carves that are essentially linear. For the cases 
In which this linearity of the lift curves Is not 
present - , for example/ near the oolnt at which noraal 
stall or premature coitpresslblll ty stall occurs - the 
direct method is no longer satisfactory. For such 
oases, a third trjethod of determining span load distribu- 
tions, referred to herein as the "generalized method," 
makes use of a orocess of successive approximations to 
apply the llftlng-llne theoi^r. Applications of this 
method are giver In references 1, 2, and 5* 

The discussion that follov/p is ooncomed with the 
atrip me bhod, the conventional method, and the generalized 
method of determining span load distribution with par-, 
tloular reference to their usos as related to the span 
load distribution that mdy be expected at high speeds. 

Strip Vethod 

In the strip theory, the assumption is that the 
wing Is made up of airfoil strips, each of which maintains 
its two-dlmjenslonal ^ or Inflnlte-aspeot-ratlo) lift 
characteristics. For a selected wing angle of attack 
the lift coefficient at each span station Is determined 
from the section data by ploking off the local lift 
coefficient corresoondlng to tne geomotrlo angle of 
attack at that section. Inasmuch as infinite aspect 
ratio Is assumed, the effect of varying spanwlse dis- 
tribution of Induced angle of attack In neutralizing 
any discontinuities and high gradients in the loadings 
Is neglected. The loadings are derived by multiplying 
the .lift coefficients at each station by the ratio of 
the chord at the station to the mean chord; thus, a 
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load coefficient 07— la defined, which I3 related to 

'c 

the usual wing lift coefficient G-^ as follows: 

The application of the strip theory requires no further 
explanation as it is one of the .oldest methods used in 
predicting load distribution. 

Conventional Method 

The conventional method of applying llftlng-llne 
theo3?y takes Into consideration the effect of aspect 
ratio (the -presence of downwash) ; it Is made adaptable 
for computing purposes by the assumption that the In- 
dividual span stations have constant lift-curve slopes. 
Inasmuch as airplane designers are in general familiar 
with both the Fourier series method and the Schrenk 
approximate mathod as given and discussed in reference Ij., 
no details of this method are given herein. 

Generalized Tifethod 

The need for nore generalized methods of applying 
llftlng-llne uheory (references 1, 2, and 5) arose from 
the fact that methods were desired that onuld treat 
cases In T;hi3h either partial stall or nonllnearity In 
the lift curves resulted. Tn principle, these methods 
are straightforward; that Is, fron-i the fundamental down- 
wash equatlon-a spanvlse distribution of downwash angle 
Is found for some initial assumed loading and, from the 
differences between the goometrlc and the computed down- 
wash angles at each station of the span, the effective 
angles of attack are determined. When the effective 
angles of attack are applied to each section lift curve, 
lift coefficients at each station are obtained which, 
when multiplied by the ratio of the chord at the station 
to the mean chord, define a new "check" distribution. 
The second assiimed soan loading may be taken between 
the first approximation and the check points (not neces- 
sarily a mean). The process Is continued until the cbeck 
loading coincides with that from which it was derived. 

Because of the difficulty In evaluating the down- 
wash angles, the methods of successive approximations are 
generally tedious. The downwash angle at a span station 
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Ib obtained by. carrying through the operations Indicated 
by the fundamental equation for the downwash angle at a 
spanwlse point, namely 



The pT'lme la used to indicate the span point at which the 
downwash angle Is belnQ found. 

Reference 1 carries out the Indicated ooeratlon by 
plotting the load coefficient 3^|- against l/(y' - y) 

and Integrating the area to obtain the downwash. A 
taofrS' coirplloated method given In reference 2 carries 



out the Indicated operation by olottlng — ^ — 

y - y 

agalnsb the span station y, and Integrating to find the 
area.' In both methods, since the denominator becomes 
zero when y = yS the Integration mast exclude a small 
area on either side of the singularity, which Is then 
soperately oonsldered by en approximate forin'.ila. Of 
these two methods, that of reference 1 Is simpler bocause 
It does not depend upon the grcphlcal determination of 
the. slope of the loading; curve. A much less laborious 
method was developed In referonce 3* which Is demonstrated 
In reference 5* ^^^^ has been fo\ind by experience to entail 
only a small fraction of the labor of the methods of ref- 
erences 1 and ^. 

The method of reference 5 Is essentially as follows; 
The distribution of span loading Is expressed^ by the 
general formula (from reference 6, equatl-jn (iil)) 




• (1) 





(2) 



r=0 



7 



which Is then substituted Into the basic downvash- angle 
fozmila (equation (1)) and Integrated step by step to 

. . arrive at the series 

■ ' — -- - , 

f ' = 2- ^2r 82r (?) 

r=0" 



where 
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2r 




2n + 1 



P - rP 
*^r-s "^r 



-s-lj 



r2S 



and 

P = 1x3x5 (2r - 1) 
^ 2xii.x6 . . , 2r 



At this point In the development glv&n In reference 5 It 
Is assumed that this downwash angle iray be expressed as 

and since, by equation (2), Oj^^ may be given In terms 

c 

of j^, yj- • • y^' 

^ = Yi - ""2' 

m=l r=0 

When ■ the assumed downwash angle given by equation ( 5) 
Is onmpared wltli the downwash angle given by equation (3)* 
the corsespondlng terms for the two series may bq made 
to agree If y^, yg, • • • ^ \» ^2' * ' * ^n 
satisfy the condition 

ym^"" - ym^"* = 82^.^^ ( 6) 

m=l 
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Inasmuch as equation (3) contains only y^^, 

^2' • • • ^1* ^2' • • • ^n' values of y 

may be prescribed Independently of the a^r coef- 

floienta and, therefore , Independently of the loading 
^0^,=^; the problem consequently reduces to the solu- 
tion of a system of simultaneous equations. 

The constants for the case of five siiriultaneoua 
equations (n = 5) as taken from reference 3 are given 
in the following table and have been determined to 
yield the downwash angles in degrees: 



X 


0 


0.5 


0.5 


0.7 


0.9 




150.10 


-72.67 


-16.16 


21. 1;2.: 


-107.60 




-1214..14.6 


187.77 


-75.50 


-Si. 12 


2ij2.99 




8. 85 


-78.17 


200.1k 


-io.i|.3 


-211.G.86 


\ 


-l2.i^9 


-10.52 


-76.55 


169.98 


38.81 


^5 




-5.96 


-8.qo 


-63.6l^ 


211.65 



Tf equation (Jj) is expanded to n = 5» downwash- 
engle equation becomes 




9 



and, for example, for y = O.J 

Early experience gained In applying the methods of 
references 1, 2, and 3 has indicated that Identical 
results would be obtained with each of the luethods 
after a sufficient number of trials but the method of 
reference J ''^3 found to be much easier than the others 
and consequently was the method used herein. A tabular 
scheme was finally worked out, which enabled the span 
load distribution to be obtained very quickly. A 
sample of the computations Is shown In table I. 



DISCRIPTION OP t.T:NGS AtTD BASIC DATA 



The hypothetical wings (see fig. 1) foi- which the 
high-speed loadings were determined were both of ellp- 
tlcal plan form and had a thlcknGss ratio varying from 
17 percent of the chord at the root to 9 percent at the 
tip. One wing was assumed to be of NAG A OOXX sections 
and the other of hlgh-crltlcal-speed 16-5XX sections. 
Tbe notatJon XX Is substituted for the airfoil thick- 
ness In the designation. A linear geometric twist that 
would most nearly Introduce zero aerodynamic twist at a 
Mach number of O.JO was applied to the wings. In order 
to accomplish this result the 16-5XX wing was given a 1° 
geometric washout, whereas the symmetrical OOXX wing 
fulfilled the condition with zero geometric twist. 

The basic data for the OOXX airfoil sections were 
obtained from the British National Physics Laboratory 
and those foy* the I6-5XX sections from the NAG A 2J4.-inch 
high-speed tunnel. The pertinent facts conoemlng the 
scope of the data and the reference papers in which the 
material is presented are given in the following table: 



BASIC DATA 



Airfoil 
section 


Angle-of- 
attack 
range 
(deg) 


Chord, 
(in.) 


Kach number 
range 


Approximate Reynolds 
number range 


Ref- 
erence 


NACA OQXX sections 


0012-63 
0015-63 
0017-63 
0020-63 


"h to 12 
0 to 5.75 

0 to 5'5 
-2 to lli. 


2 

I.U95 

1.5 
1.2 


oJiO to 0.725 
.30 to .80 
.30 to .60 
.30 to .60 


k50,ooo to 7Go,ooo 

350,000 to 500,000 
360,000 to 500,000 
200,000 to 500,000 


I 

9 

10 


NACA 16-5XX Elections 


16-506 

16-509 
16-512 
16-515 
16-521 


-3.2 to 2 
-3.6 to 3 
-k to I|. 
-2.3 to L..5 
-1.6 to h 


c 

y 

5 
5 

c 
J 

5 


0.30 to 0.75 
.30 to .75 
.50 to .75 
.30 to .75 
.50 to .70 


700,000 t6w2*Q00,000 
700,000 to 2,000,000 
700,000 to 2,000,000 
7X,000 to 2,000,000 
700,000 to 2,000,000 


11 
11 
11 
11 
11 



The baslo Beotiozi data were so pl-dtted as to permit 
nonlinear interpolations of th^,JL±^t curves for any of 
the airfoil 'seotions along, the wing. ..The unusual type 
of plot adopted vas suggested hy one used in' reference 10 
for a different purpose and was constructed as follovs; 
For a given Mach number the usual section lift curves for 
various thickness ratios were plotted on the same figure. 
The angle-of-attack scale for each thickness was staggered 
a distance protiortional to the airfoil thickness (heavy 
lines in figs, 'd to 6) and the points of equal angle of 
attack were then Joined, producing a surface or "carpet" 
effect. This surface allows the lift curves to be. 
easily Interpolated for airfoil thiclcnesses between 
those represented by the basic data and, in addition, 
gives a pictorial representation of the airfoil family 
making up the wing. 

The particular airfoil seotions used were chosen 
because the data were available from the same wind tunnel 
for each family. The 16-5XX sections were tested through 
the sdme Reynolds number range; however, the COXX air- 
foils, as shown in the tabulation of the data, were of 
different chords and consequently the Reynolds numbers 
were not the same. Also, according to reference 8, the 
OOl^-S^ airfoil section was obtained by cutting down 
the 0017-63 section and the erd blocks were probably 
twisted ona-fourtli of 1° relative to the new datum, 
giving an observed error in fhs results in the aero- 
11ft anfrle. Nevertheless, no attempt was made to cor- 
rect the basic data since it was considered that the 
more or less arbitrary corrections v;ould not alter the 
results an amoua.t sufficient to warrant tsmoerlng with 
the basic data. 

As shown in the tabulation of the basic data, the 
Mach number range for the OOI2-63 airfoil section was 
0,kO to 0.725. The data for this section at M = O.3O 
(fig. 2) are therefore shown as extrapolations. In 
the carpet for the case of M = O.75 (fig. 6) the data 
for the 12-peroent- thick section, however, are shown as 
neither basic nor extrapolated data but rather as a 
heavy dashed line to indicate that the line, although 
not from basic data, was better than the extrapiol ated 
curves because the extension of the data from O.723 to 
0.75 slight and was made by meanB of an auxiliary 
carpet. In this carpet, which is not shown, the lift 
curves were plotted for the various Mach numbers of the 
0012-63 airfoil and the- curves were then spaped In pro- 
portion to the Mach niunber increments. 
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RESTJLTS OP COMPUTATIONS 



Span load distribution computations using the' strip- 
theory inethod end the generalleed method of applying 
llftlng-llne theory were perfosrmed for each of the hypo- 
thetical wings at five or six' angles of attack of the 
root section and for each of the five Wach numbers for 
which data were available, that Is, from figures 2 to 6. 
The hypothetical wings were purposely chosen ^Jlth ellip- 
tical plan form and with zero aerodynamic twist In order 
that the span loading which would be obtained by the 
conventional application of llftlng-llne theory (refer- 
ence Li.) woxjld be elliptical at the various lift coef- 
ficients . 

A sample of the computations required and the 
procedure followed at each value of angle of attack and 
Mach number to obtain the span load distribution by 
either the strip method or the generalized method of 
reference 5 shown In table l.t The computations 
shown epoly to the I6-5XX wing for an angle of attack 
of 2° at the root and for a yach nTimber of O.75. Tn 
this table the compute.tlons required for the strip- 
theory calculations are included in the first six 
columris of the first block. The remaining columns 
and blocks ere required to corcplete' the conputatlons 
for the method of successive approximations. It will 
be noted both in the table and In -figure 7« which 
supplenisnts the table, that In the computations by the 
■method of successive approxJiLStions the first aoproxl- 
matlon of the loading Is assumed to be that defined 
by the strip theory for which the local lift coefficients 
are shown by th^ solid line in figure 6. The check 
points 6btalned,'for the first approximation are shown 
both in figure 7 and on the carpet of fleure 6. Table T, 
in conjunction with figure 7 and figure o, may be 
followed to show completely the detailed procedure. 

The span load distributions obtained from each of 
the con^jutations were then integrated and the results 
cross-plotted against wing lift coefficient to 

obtain loadings as a function of f'ach number at a 
ntunber of evenly spaced wing lift coefficients equal 
to 0.1, 0.2, 0.3f O.Ut and O.5, The results of the 
cross plots are shown in figures 8 and 9. The loadings 
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given in these figures were obtained both by use of the 
llftlng-llne theory and by the atrip theory. OSie solld- 
Ilhe curves in these figures represent the reference (or 
conventional) design span loading, which has been as- 
sumed by designers not to change with TAach number.- 

The span loadings of figures 8 and 9 were in turn 
integrated to give a measure of the bending moment at 
the wing root end at the ^^-percent span station. The 
results of these integrations are shown in figures 10 
and 11, in which a bendlng-moment coefficient is 

plotted against the wing lift coefficient. The bendlng- 
moment coefficient shown in these figures is defined 
mathematically by the expression 

^1.0 



/. t.t) 



This definition requires that, at the wing root, the 
bendins-moment coefficient equal the vlng lift coef- 
floienb multiplied by the lateral center of pressure of 
the sp£n loed. Also, with the bendlng-moment coef- 
flclcn!: defined in thj s manner, the magnlturie of the 
bendlr]^ moment at any station y* is given by 



Bending Tfoment = — p- qSb 



DISCUSSION 



With respect to the generalized method used in the 
span load computations it may be seen that when a sys- 
tematic tabulation system is adopted, very little work 
would be required even in a general case. The table as 
given could be shortened still further because the first 
six columns in the second, third, and fourth approxi- 
mations could be deleted. All the columns have, how- 
ever, been Included in this paper in order to show each 
step clearly. Those familiar with span load confuta- 
tions will easily recognize that the amount of work is 
less than that required by the Fourier Series method 
of reference h and is only slightly more than that re- 
quired by the optional but shorter sohrenk approxi™ 
mate method that is included therein. 



.It will be noted that, althougki far more coniplete 
basic aeotlon data were available for the families chosen 
than for bther airfoil families, the data are still 
rather limited In scope both as to the Mach number and 
the llft-ooefflolent range obtained. For this reason 
It Is Impossible to Illustrate either the effects of 
yach number on span load changes near values of M = O.85 
that some modem airplanes have reached or to Illustrate 
changes that occur with even moderately high lift coef- 
ficients. Examination of the results given In any one 
of figures 8 through 11 Indicates that changes In loading 
occur rapidly for the I6-5XX wing at a Mach number above 
0.7 and for the OOXX wing at a Wach number above 0.6 and 
a lift coefficient above O.I|. The highest value of 
Mach number for which data are available is only 0.75" 
In the case of the OOXX wing, sufficient section data are 
available at this Wach mmiber to allow computations to 
be made to a wing lift coefficient of only 0.2 and for 
the 16-5XX wing to a wing lift coefficient of 0.5. 

From qualitative considerations It would be ex- 
pected that the soon load center for a wing composed of 
a ooncistent family of sections mlglit shift either in- 
board or outboard, depending upon the combination of 
Mach niimber and angle of attack at which the wing is 
oper&ting. At high ''-aoh numbejrs and low to moderate 
lift coefficients, as shown by 'ilie rtsults of the 
present analysis, the early stLll of the thicker in- 
board sections of the wing produces an outboard shift 
of the center of lorad. At relatively low Vaoh numbers 
and high angles of attack, however, the loading may 
shift Inboard because of the earlier stall of the ".■ 
sharp-nosed outboard sections. The span load center 
could also move inboard when the wing operates at low 
angles of attack and near critical I.Tach numbers, because 
the higher negative pressures over the thicker Inboard 
sections would then be expanded more than the lower 
pressures of the thinner tip sections. The basic data 
are not of sufficient range to show the Inboard shift 
of the center of load caused by tip stalling; howaver, 
the Inboard shift as caused by reduction of pressure 
over the upper surface of the root sections is indicated 
by the results given in figure 8, in which the loading 
may be noted to move inboard for the 16-5XX wing at 
lift coefficients below 0.20 at M = O.7O. 
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Fop a given rigid wing the shift In span loading 
that occurs at high Hach nvmibers depends on the change 
with Hach number of the relative positions of the lift 
otirvea for- thfr' fleotlonr< making the wing* -The.. re- 
lative positions of the section lift curves at various 
Hach numbers may be clearly noted from the carpet plots 
(figs* S to 6) ; these plots are particularly advantageous 
for presenting a quick impression of the high-speed per- 
formance of the wing. ' For Instance, for the OOXX wing 
at a Mach number of 0*70 (fig* 5) the relative positions 
of the section lift ourvew for the 9-percent and 17- 
percent-tblck sections at low angles of attack eve almost 
the same as for the ifrlng at a Mach nianber of 0*30 (fig. 2) 
and little change In loading Is therefore to be expected 
at these angles* The relative positions of the section 
lift curves for the Ixiboard and outboard sections change, 
however, as the angles of attack Increase, and the span 
loading changes should become more severe with Increasing 
angle of attack. This conclusion Is borne out In the 
bendlng-moment-coeff Iclent curves of figure 11, which 
show that the bendJ ng-moment coefficient for M = 0,70 
Increases with lift coefficient. For the 16-5XX wing 
at a Mach number of 0,75 (fig* 6), all the section lift 
curves of the thicker sections have been displaced down- 
ward. Indicating that the Inboard sections carry less 
lift but. Inasmuch as the slope changes little between 
the outboard and Inboard sections throughout the angle- 
of -attack range, the Increase In b ending-moment cioef- 
f Iclent should be constant with angle of attack. Fig- 
ure 10 verifies this deduction^ 

From the foregoln/: discussion It Is evident that the 
manner In which changes In the span loading occur Is de- 
pendent upon the airfoil sections and other geometric 
characteristics of the wing. For the two wings considered* 
the changes of span loading with Mach number brought 
about in the first case a constant bending-moment-coeff Iclent 
Increase over the low-speed value with Increasing lift 
coefficient and In the second case a gradual Increase 
oyer the low-speed value with Increasing lift coefficient. 
F-rom an examination of the carpets of different wings at 
both low and high Hach numbers, the family of airfoils may 
be selected that has the most favorable span-load- change 
characteristics* 

When related to possible flight conditions, the 
b ending-moment Increases Indicated In figures 10 and 11 
are sufficiently high that the safe load factor would be 
considerably diminished* Figures 12 and 13 have been 
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prepared to Illustrate this point. Figure 12 shows the 
variaticfn of bending-moment ooeffloient with load factor 
at two spanwlse stations for a pursuit airplane equipped 
with a 10-53DC wing having a loading of I4D pounds per 
square foot and a }.?aoh number of O.75 at 10,000 feet. 
Figure 13 shows the some results for the airplane eq\iipped 
with 'the OOXX wing at a Mach number of O.70 at 10,000 f^et. 
Under 6rdlnary conditions the wings would be designed to 
sustain the bending mo" ents corresponding to a load 
factor of- 8g and bendiiig-moment coefficients of 0.0i|.5 
and 0.-255 the percent and root stations, respec- 
tively. Both the strip theory and the generalized 
'lifting- line theory show that these valvies of bending- 
moment coefficient would be obtained at values of the 
load factor substantially less than 6g. 

In the strip theory, no consideration is given the 
effect oif induced flow in leveling the load gradients 
along thd span; the loads obtained with this method are 
therefore more severe than those obtained when the 
generalized lifting-line theory is used. For the two 
hypothetical wings, however, even this simple strip 
theory gives results closer to the high-speed loading 
and is more conservative than is the conventional method 
used In design. 

- It should be noted that the downwash equation has 
been assumed to apply for compressible as well as for ■ 
incompressible flow. To a first approximation this 
assumption is reasonable because, regardless of the type 
of flow in which a given lift is realized, the principle 
of Induction would -still apply. Even substantial 
changes in the downwash angle due to compressibility 
should', however, not be- expected to greatly alter the 
span loadings, as these changes would bccur all along 
the span and, for- a given family of section lift curves, 
it Is the" change of .the form of the downwash-angle 
distribution along the' span which Influences the span 
load distribution moat. • 

The results given heroin were obtained for a rigid, 
smooth wing, without fuselage and nacelles. In an 
actual cAse the crit?.cal condition may occur at lower 
yach numbers because of flow disturbances caused by 
protuberances near the root station, such as inspection 
plates 9nd wing fold doors, and by fuselages and 
nacelles. 



17 



.The effect of i?lns naoelle and fuselage Interferenoe 
on the span loading may be taken into account , to a first 
order at least, by considering that these items cause a 
, - .change in .the .Macli.nujnjbej^a., Qf..the sections along the wing 
span. This effect dioilnlshes as the distance is in- 
creased from the interfering body. Reference 12 indi- 
cates a method by which the increments in Maoh number 
caused by interfering bodies such as nacelles and fuse- 
lages may be determined. The application of computed 
increments in a practical case, however, suggests that 
the basic data given in the carpets of figures 2 to 6 be 
plotted to obtain five- charts, one for each thickness 
at the selected spanwise .stations, with the scales off- 
set to give the Mach number instead of the thickness as 
the parameter on each phart. The lacTf and uncertainty 
of section data at hi^h soeeds on various families, how- 
ever, make the carrying out of such a detailed procedure 
hardly warth while at the present time. 



CONCLUDINa REMARKS 



The spanwise dist- Ibution of load on an airplane 
wing at high speeds may be determined by means of a 
generalized method of applying llftlng-llne theory. 
The results of applying such a method show that the 
bending-moment changes that can occur at high Vach 
numbers may be sufficiently great to render necessary 
a modification in the procedure now used in computing 
span loading. 

In order to determine the validity of the general 
application of the lifting-line theory, it is recommended 
that measurements of the span load distribution be made 
in a high-speed wind tunnel on a wing coii?iosed of airfoil 
sections for which the two-dimensional high-speed char- 
acteristics ere available at fairly high Mach numbers and 
at high angles of attack. 



Langley Memorial Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Langley Field, Va, 
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Figure 5. - Carpets ^r /nfcrpo/af/on of bos/c section data of f6-6X/ancta)X)C 

airfo/ts. tW=070 
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F/gure 6 .- Carpets for mtcrpc^fion cf das/c S€cfyor)clata /6-5XX and COXX 

atffo//5. M=0.73. 
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(a) Lifting-linc fhcory y0^i9r>_ Cd) Strip theory 
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Figure 6 rTfic effect of A/facfi number or? span foaof/ng for 
various lift coefficients on on ffACA I6-5XX 
Wing. 
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Figure 9 r Tfie effect of Mach r)umber on span toad in g 
for vanous f/ft coeff/c/cofs on an NACA OCX/ 
winq. 
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